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Abstract 

The dependence of the interwall conductance on distance between walls and relative positions of walls are calculated at the low voltage by 
Bardeen method for (n, n)@(2n, 2n) double-walled carbon nanotubes (DWCNTs) with n = 5, 6, . . . , 10. The calculations show that interwall 
conductance does not depend on temperature (for T ^ 500 K) and current- voltage characteristic is linear The conductance decreases by 6 orders 
of magnitude when the interwall distance is doubled. Thus, depending on the interwall distance, DWCNTs can be used as temperature stable 
nanoresistors or nanocapacitors. 



1. Introduction 

Carbon nanotubes (CNTs) have attracted much attention from the viewpoints of basic science and technological applications 
because of their unique electronic and mechanical properties [1]. A set of nanoelectromechanical systems (NEMS) based on relative 
motion of carbon nanotube walls in the multi-walled CNTs have been proposed (see Ref. [2] for a review). Some of these NEMS, 
such as variable nanoresistors [3,4], are based on tunneling current between CNT walls. The operation of these nanoresistors have 
been considered on example of double-walled carbon nanotubes (DWCNTs). DWCNTs produced by high-temperature treatment 
of single-walled nanotubes with fullerenes inside [5,6], chemical vapor deposition [7,8] or arc discharge in hydrogen [9,10] have 
interwall distances ranging from 0.33 to 0.42 nm. However the previous calculations of tunneling current between DWCNT walls 
were done only for DWCNTs with interwall distance equal to 0.34 nm [4,11]. 

The another type of nanotube-based NEMS is a nanocapacitor with coaxial nanotube walls used as its plates (for example, 
nanotube-based variable nanocapacitor can be used for controlling of operation of three-terminal memory cell by electrostatic 
forces [12]). In such a nanocapacitor the interwall tunneling current is a leakage current and interwall distance should be chosen so 
that tunneling current is negligible. Thus calculations of interwall tunneling current as a function of interwall distance is an actual 
problem for elaboration of nanotube-based NEMS. 

In this Letter, we calculate the dependence of the interwall conductance in DWCNTs with armchair walls on distance between 
nanotube walls at the low voltage using the Bardeen formalism within non-orthogonal tight-binding approximation. The influence 
of temperature on the interwall conductance is investigated. According to semiempirical calculations [13-15] and calculations based 
on density functional theory [16] the interwall interaction energy of DWCNTs with armchair walls depends on the relative positions 
of walls. Here we show that the total transmission and therefore conductance of DWCNTs with armchair walls depends also on the 
relative positions of walls. 

2. Tunneling matrix element in DWCNTs 

In the tight-binding approximation, the wave function of graphene can be expressed as [17] 

where fc is a two-dimensional vector in reciprocal space of graphene lattice with lattice constant a ~ \/3ac-c (ac-c is the C-C 
bond length), d = (oi + a2)/3 is the vector between two atoms in the unit cell, ijj{k) = 1 + exp(— ifcai) + exp(— iA;a2), 
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Fig. 1. Geometry of the (n, n)@{2n, 2n) DWCNT corresponding to tlie relative position of the walls with A^l? = and AZ = 0. Primitive unit cell (shaded area) 
is represented in the lower part in which the dashed line con'esponds to the inner wall projection on the external one (solid line). The number of coincidental atoms 
in the primitive unit cell is 2 and in the DWCNT translational unit cell is = 2n. The value = 2n repeats with periods ip^ = it/10 and = a/2 at the 
relative rotation of walls and their relative displacement along the DWCNT axis (in z direction), respectively. 

Rg ~ gidi + here ai and 02 are the lattice vectors of graphene, and Nq is the number of the graphene unit cells. The 

function is the Slater 2px orbital [18]. The signs ± correspond to the tt and tt* configurations of orbitals in graphene. 

For (n, n)@(2n, 2n) DWCNTs rotational period of the inner wall 27r/n equals to 2 rotational periods of the external wall, and 
their translational periods are coincide and equal a. Therefore the structure of the DWCNT is determined by the structure of the 
primitive unit cell restricted by cylindrical coordinates ^ 1^9 < 27r/n, ^ z < a. These DWCNTs are commensurate and belong 
to the "strong" coupling case in the terminology of [19]. 

An ID Bloch function of (n, n) CNT in state m for a 2px orbital xa,b {A and B in Fig. 1 correspond to two nonequivalent atoms 
of graphene unit cell) in cylindrical coordinates (p, (f, z) is (see Ref. [20]) 

P> z) = —= ^ exp(iA:Zj) ^ exp(im$;)XA,B(p - Pw, <^ - ^hZ - Zj), 



^ j=i 1=1 

where k is the wavenumber of an electron along CNT, m = 0, 1, . . . , n — 1 is the subband number, Zj and — 2TTl/n are 
cylindrical coordinates of carbon atoms, is the radius of the CNT, and N is the number of nanotube translational unit cells. The 
second sum corresponds to rotational symmetry of the CNT. 

In the Bardeen formalism [21], the tunneling matrix element between states \I/in and \I'ex of the inner (in) and the external (ex) 
walls of a DWCNT, respectively, is 

Afi„,ex = / (*rnV*ex - *exV** ) • dS , 

2too Js 

where S is an arbitrary cylindrical surface with radius ps (see Fig. 1) between the inner and the external walls, h is the Planck 
constant, and toq is a mass of the electron in vacuum. 

The matrix element of the probability density current operator J between states ^Pin and \I>ex of the inner and the external walls 
for Xin and Xex orbitals is 

1 



in, ex 

N N 



exp[i{ki„Zj - fcex^j')] ^ ^exp[i(TOin^'; - mex^*;')] 



X (Xex(P - Pex^ip- Zj,)\J\Xm{p " Pm,f- ^l, Z - Z j)) , 

where pi„ and pex (see Fig. 1) are the radii of the inner (in) and the external (ex) walls, respectively. 

Figure 1 shows that for (n, n) @ {2n, 2n) DWCNTs in the case where one or two atoms of the primitive unit cell of the inner wall 
have the same cylindrical coordinates with some atoms of the outer wall = and Zji = Zj, the differences |$/' — <i>;| and 
I Zj/ — Zj I for other pairs of atoms are essential and their overlaps are small. Therefore only the interaction between 2px orbitals of 
atoms with the same cylindrical coordinates $/' = $/ and Zji = Zj is considered. In this case the matrix element can be rewritten 
as 



M:2T'^ = J2 exp[i(fcin - fcex)Zj)] J2 exp[i("im - mex)$i)] 

"'"^^ J=l 1=1 

X (Xex(p - Pex, (p-<^l,Z- Zj)\J\xin{p - Pin, (p-<^l,Z~ Zj)) 



2 
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Fig. 2. Configuration of interaction between two states n and n* for atoms witli the same cylindrical coordinates $ and Z in DWCNT. 

= 4,„,.„<^™,„.™„^ (X* ^Xex - Xex^X* ) dS, (1) 

where (5fc.„,fc„ and (5m.„.„i„ are Kronecker symbols ensuring impulse and z projection of orbital momentum conservations in the 
electron tunneling process, respectively. 

We take p5 — (pin + Pex)/2. In this case (dxex/dp)p=ps = — (dxex/d/o)p=ps, and the matrix element is 



™o Js 



(2) 



where X = (CV327r)^/^ ^ exp(-^r/2) is the Slater orbital; ^ = 1.5679/aB [18], with Ob is the Bohr radius, and r is module 

of the radius vector from atom centre. 

The configurations of interaction between two states tt and tt* for the atoms with the same cylindrical coordinates $ and Z in 
DWCNT are shown in Fig. 2. In the case of tttt* interaction (Fig. 2(b)) the matrix element is zero because of opposite signs for A 
and B interacting orbitals. 



3. Interwall conductance in double-walled armchair carbon nanotubes 



When a voltage U is applied between the inner (in) and the external (ex) walls, the tunneling current / can be written as [22] 

n 



in ex 

where e is the electron charge module, and /(£'in) and f{Ee^) are the Fermi-Dirac functions. The delta function 6{Ein — E'ex) 
ensures the energy conservation for the tunneling between walls. The summation is carried out over k and m quantum numbers of 
walls. 

The energy dispersion relations E{k, m) for the armchair {n, n) CNTs are written as [1] 



/ / imT\ / ka\ „ f ka\ 

E{k,m) = ±7Jl±4cos( — j cos( y ) + 40082(^1, (4) 

where 7 — 2.6 eV is resonance integral [23,24], k is the wavenumber of an electron along the CNT, and m = 0, 1, . . . , n — 1 
is the subband number. According to Eq. (4) for (2n, 2n) CNT the distance between Ep ~ and minimum of m = ±1 band 
is eUo = 7^/1 — cos2(7r/2n) w ^n/2n. The band structure of the (5,5)@(10,10) DWCNT is presented in Fig. 3. In this case, 
eUo w 0.8 eV and for bias \U\ < Uq for temperature T ^ the tunneling occurs only between states with min = rriex ~ 0. In 
this voltage domain, the energy bands E{k, 0) are linear functions of k and density of states is constant. Direct calculation of the 
fc-integral in Eq. (3) for the armchair nanotube of length L ~ Na gives (see Appendix) 

/ = ^|i^CJ2^,„(£;p)Ax(^;F)i^ (5) 

where M^^^^ is the tunneling matrix element for the translational unit cell in the case mi,, = mex = 0, Di„{Ep) and D^^{E-p) are the 
densities of states per unit length along the nanotube axis at the Fermi level E^ — for the inner and the external walls, respectively. 
The formula (5) shows that the tunneling current at \U\ < Uq does not depend on temperature when k^T < O.leC/o (O.leUo for 
n = 5, 6, . . . , 10 corresponds to 500 K) and current-voltage characteristic is linear 

Formula (5) is also valid for any non-constant £'in(ii'F) and Dg^{Ep) in the low temperature limit (T — > 0) and small voltage 
(compare with Ref. [22]). For armchair Din{Ep) = Ds^{Ef) = 8/(7170-^8), i.e., the interwall tunneling conductance per one 
translational unit cell (L ~ a) is 
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Fig. 3. Schematic band structure of the (5,5)@(1(). 10) DWCNT in the equilibrium state (left) and under bias \ U\ < Uq (enlarged); fcp = 2n/3a and Ep = 0. 

We have studied the influence of the relative position of walls of the DWCNT on the tunneling matrix element. This position is 
determined by the angle A<i> of the relative rotation of walls and the relative displacement AZ of walls along the nanotube axis, 
see Fig. 1. For a discrete set of relative positions of walls (A$c, AZ(,) some atoms of the inner and the external walls have the 
same cylindrical coordinates $ and Z. The number of these "coincidental" atoms in translational unit cell of (n, n)@(2n, 2n) 
DWCNT (which consists of n primitive unit cells) can be written as 

( 2n, if (A$(;, AZc) ~ {qn /lQ,pa/2), where q andp are integers, 
\ n, if (A<i>(;, AZc) ~ ((77r/30,pa/2), where q is not multiple of 3. 

For relative positions of walls (AZ, A(l>) = (AZc, A^c), we numerically calculate the transfer matrix element for one pair of 
orbitals ^™ ^ = M°^,Jnc using Eq. (2). The calculations for the (5,5)@(10,10) DWCNT give J = 0.325 eV for one pair 
of Px orbitals. This value is approximately equal to the value used in Ref. [20]. 

For small relative displacements of walls from the relative positions (A$, AZ) = (A$c, AZc) we neglect the "broadening" of 
J-s in Eq. (1) and calculate dependence of the matrix element of the (5,5)@(10,10) DWCNT on A$ and AZ according to Eq. (2) 
taking into account the dependence of x on ™d AZ. Figure 4 shows the dependences of square of the tunneling matrix element 
per one translational unit cell and the conductance G measured in quantum units /ttK at 7 = 2.6 eV on A^/ip^ and 

AZ/flc, where — tt/IO and — a/2 are the angular and translational periods of these dependences for the (5,5)@(10,10) 
DWCNT, respectively (Fig. 1). 

Note, that the number of coincidental atoms (for A$ = 0, AZ — 0) per one translational unit cell in (n, n) @ (2n, 2n) DWCNTs 
is 2n, which must be taken into account in calculating the tunneling current in this case. Figure 5 shows the calculated dependence 
of the interwall current I /Id for A$ = and AZ = on the interwall distance pex — Pin at the low applied voltage, where / and 
/(J are the interwall currents in [n, n)@{2n, 2n) and the (5,5)@(10,10) DWCNTs, respectively, with n — 5, 6, . . . , 10. It is seen 
that the ratio of the interwall currents of DWCNTs with armchair walls decreases nearly exponentially with increase of the distance 
between the nanotube walls. The line in Fig. 5 corresponds to the fit log(///d) ~ 15.4 — 1.94(/9ex — Pin), where pex and pia are 
measured in nanometers. 




Fig. 5. Dependence of ratio of the interwall currents ///j on the distance 
Fig. 4. Dependences of square of the tunneling matrix element \M^^^^^'\^ and between walls (pcx — Pin), where / and are the interwall currents of the 
conductance G for the (5,5)@(10,10) DWCNT per translational unit cell on (n, n)@(2n, 2n) and (5,5)@(10,10) DWCNTs, respectively, pi„ and pcx are 
the relative position of walls, where ip^ = tt/IO and = a/2 are the angular the radii of the inner and the external walls in a DWCNT. The straight line is 
and translational periods of these dependences for the (5,5)@(10, 10) DWCNT. the fit of values. The calculated results correspond to A<I> = and AZ = 0. 



4 



r " 



DWCNT 



(n, n) CNT 



(2n, 2n) CNT 




Fig. 6. Scheme of a DWCNT-based NEMS which (depending on the interwall distance) can be used as a variable nanoresistor or variable nanocapacitor 

4. Conclusion 

In conclusion, the interwall conductance in the DWCNTs has been studied using the Bardeen method. Numerical calculations 
show that in the (5,5)@(10,10) DWCNT relative sliding of the walls along the DWCNT axis leads to 30% variation of the interwall 
tunneling current, while the dependence on the angle of relative wall rotation is more pronounced. It is found that conductance 
of {n, n)@{2n, 2n) DWCNT for n = 5, 6, . . . , 10 does not depend on temperature for T ^ 500 K. The interwall conductance 
G of (n, n)@{2n, 2n) DWCNTs decreases by 6 orders of magnitude when the interwall distance is increased from 0.34 nm to 
0.68 nm. Thus DWCNTs can be used as temperature stable (up to 500 K) nanoresistors in the case of small interwall distances and 
as nanocapacitors with negligible leakage current in the case of interwall distances ^ 1 nm. The scheme of these NEMSs is shown 
in Fig. 6. 

Acknowledgements 

This work is supported by the Belarusian Foundation for Basic Research (Grant No. F08R-061) and Russian Foundation for 
Basic Reseai-ch (Grants Nos. 08-02-90049-Bel and 08-02-00685). 

Appendix A. Calculations of the tunneling current 



In the case — m^^ — 0, from Eq. (4) we obtain Ei„{k, m = 0) = E^,^{k, m = 0) — i?o(fc) — ±7[1 — 2 cos(fea/2)], and the 
formula (3) for the tunneling current in DWCNTs with the length of L (using = (^/^tt) / dfc) can be written as 

IiU,T) _ 64e 



L2 



dfci„ / dfcex{/[£^o(fcin)] - f[Eoih,) + eU]}\^Cj'6[Eo{h„) - So(fcex)]4i„.fc„ 
^TTft ./bz Jbz 

-1 / 1 \ -1 



1) 



-^IACxP(^) / {/[^^o(U]-/[i^o(U + eC/]}dM„ = ^|ACxP(^)' A, (A. 

where A = / {/[i?o(fcin)] — f[Eo{kin) + ell]} dfcin, = 27r/3a, and Ef = 0. The integral is over the first Brillouin zone (BZ). 
Jbz 

Factor 64 stands for summing over vr and tt* bands with min = mgx = 0, two wave vector directions and two spin directions of an 
electron for both CNTs. 

In order to calculate the integral A, we assume that Eo{k) is linear in the domain of integration and use variable k' — —. — — — x 



dEoih, 



knT 



dki„ 



= {kin — kp)a, therefore, dk' = adku 

k[n — kp 

The integral A can be rewritten as 



A 



1 



+ 00 



dfc' 



1 + exp k' 



+ 00 



dfc' 



where C — exp(fc'). 
Hence 



a 



-In 



1 + C 



1 + exp{eU/kBT + k') 

, /l + Ccxp(eC//fcBT) 
In ' 



dC 



dC 



C(l + C) 7o C[l + Cexp(ef//fcBr)];' 



c ; V c 

Inserting (A.2) into (A.l), we obtain formulae (5) and (6). 



-fdEo; 

^in / k,„=kr 



(A.2) 



References 



[1] M.S. Dresselhaus, G. Dresselhaus, Ph. Avouris (Eds.), Carbon Nanotubes: Synthesis, Structure, Properties, and Applications, Springer, Berlin, 2001. 
[2] Yu.E. Lozovik, A.M. Popov, Phys. Usp. 50 (2007) 749. 



5 



[3] Yu.E. Lozovik, A.M. Popov, Fullerene Nanotubes Carbon Nanostruct. 12 (2004) 485. 
[4] Q. Yan, G. Zhou, S. Hao, J. Wu, W. Duan, Appl. Phys. Lett. 88 (2006) 173107. 
[5] B.W. Smith, D.E. Luzzi, Chem. Phys. Lett. 321 (2000) 169. 

[6] S. Bandow, M. Takizawa, K. Hirahara, M. Yudasaka, S. lijima, Chem. Phys. Lett. 337 (2001) 48. 

[7] H. Zhu, C. Xu, B. Wei, D. Wu, Carbon 40 (2002) 2023. 

[8] W. Ren, F. Li, J. Chen, S. Bai, H.M. Cheng, Chem. Phys. Lett. 359 (2002) 196. 

[9] J.L. Hutchison, N.A. Kiselev, E.P. Krinichnaya, A.V. Krestinin, R.O. Loutfy, A.P. Morawsky, V.E. Muradyan, E.D. Obraztsova, J. Sloan, S.V. Terekhov, D.N. 

Zakharov, Carbon 39 (2001) 761. 
[10] Y. Saito, T. Nakahira, S. Uemura, J. Phys. Chem. B 107 (2003) 931. 
[11] D.-H. Kim, K.J. Chang, Phys. Rev. B 66 (2002) 155402. 

[12] A.M. Popov, E. Bichoutskaia, Yu.E. Lozovik, A.S. Kulish, Phys. Status Solidi A 204 (2007) 1911. 
[13] A.N. Kolmogorov, V.H. Crespi, Phys. Rev. Lett. 85 (2000) 4727. 

[14] R. Saito, R. Matsuo, T. Kimura, G. Dresselhaus, M.S. Dresselhaus, Chem. Phys. Lett. 348 (2001) 187. 

[15] A.V. Belikov, Yu.E. Lozovik, A.G. Nikolaev, A.M. Popov, Chem. Phys. Lett. 385 (2004) 72. 

[16] E. Bichoutskaia, M.I. Heggie, A.M. Popov, Yu.E. Lozovik, Phys. Rev. B 73 (2006) 045435. 

[17] R. Bamett, E. Demler, E. Kaxiras, Phys. Rev. B 71 (2005) 035429. 

[18] E. Clementi, D.L. Raimondi, J. Chem. Phys. 38 (1963) 2686. 

[19] M.A. Tunney, N.R. Cooper, Phys. Rev. B 74 (2006) 075406. 

[20] R. Saito, G. Dresselhaus, M.S. Dresselhaus, J. Appl. Phys. 73 (1993) 494. 

[21] J. Bardeen, Phys. Rev. Lett. 6 (1961) 57. 

[22] J. Tersoff, D.R. Hamann, Phys. Rev. B 31 (1985) 805. 

[23] J.W.G. Wildoer, L.C. Venema, A.G. Rinzler, R.E. Smalley, C. Dekker, Nature 391 (1998) 59. 
[24] T.W. Odom, J.-L. Huang, Ph. Kim, Ch.M. Lieber, Nature 391 (1998) 62. 



6 



